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Abstract. We constrain the parameters of dynamical dark energy in the form of a classical 
or tachyonic scalar field with barotropic equation of state jointly with other cosmological 
parameters using the following combined datasets: the CMB power spectra from WMAP7, 
the baryon acoustic oscillations in the space distribution of galaxies from SDSS DR7, the 
power spectrum of luminous red galaxies from SDSS DR7 and the light curves of SN la 
from 2 different compilations: Union2 (SALT2 light curve fitting) and SDSS (SALT2 and 
MLCS2k2 light curve fittings). It has been found that the initial value of dark energy 
equation of state parameter is constrained very weakly by most of the data while the other 
cosmological parameters are well constrained: their likelihoods and posteriors are similar, 
their forms are close to Gaussian (or half-Gaussian) and their confidence ranges are narrow. 
The most reliable determinations of the best-fit value and la confidence range for the initial 
value of the dark energy equation of state parameter are obtained from the combined datasets 
including SN la data from the full SDSS compilation with MLCS2k2 fit light curves. In all 
cases the best-fit value of this parameter is lower than the value of corresponding parameter 
for current epoch. Such dark energy loses its repulsive properties and in future the expansion 
of the Universe changes into contraction. 

We also perform a forecast for the Planck mock data and show that they narrow signifi- 
cantly the confidence ranges of cosmological parameters values, moreover, their combination 
with SN SDSS compilation with MLCS2k2 light curve fitting may exclude the fields with 
initial equation of state parameter > —0.1 at 2a confidence level. 
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1 Introduction 

The discovery of the accelerated expansion of the Universe has led to the introduction of a 
new mysterious component - dark energy. Its unknown nature is one of the main puzzles 
of modern cosmology. Since the simplest explanation - a cosmological constant, faces many 
interpretational problems, a variety of alternative models have been proposed (see reviews [1]- 
[8] and books [9]- [10]). The simplest alternative approach treats dark energy as a scalar field 
with standard Lagrangian. The model is defined by the scalar field potential which can be 
either physically motivated or obtained via reverse engineering from the variables describing 
dark energy in a phenomenological fluid approach: the energy density and equation of state 
(EoS) parameter w. The latter can either be constant or vary in time. The character of 
the temporal variation of w is usually assumed ad hoc. Nevertheless, physically motivated 
dependences of the equation of state on time are sought. 

The simplest and most widely used Lagrangians of a scalar field are the Klein-Gordon 
(called also classical) and the Dirac-Born-Infeld (often called tachyon) ones. 

In our previous papers [11]- [13] we have analyzed the parametrization of the EoS by its 
current wq and constant adiabatic sound speed c^, which corresponds to the EoS parameter 
at the beginning of expansion We- In this case the dark energy EoS is of the generalized 
linear barotropic form. Such a parametrization is easy to motivate physically. For the dark 
energy in the form of the scalar fields with barotropic EoS the analytical solutions for the field 
variables and potentials exist for both classical and tachyonic Lagrangians. In the case of 
these dark energy models the relation between the current and early values of EoS parameter 
determines two drastically different scenarios for the future evolution of the Universe. 

The only way to verify the plausibility of a dark energy model is to confront its pre- 
dictions with the observational data and to find the allowed ranges of its parameters (for 
discussion of the cosmological parameter estimation see e. g. [9]). For this purpose the 
Markov Chain Monte-Carlo approach is widely used. In the paper [11] we have found that 
the early EoS parameter We remains unconstrained by two combined datasets including the 
recent data on CMB anisotropy, large-scale structure of the Universe and light curves of 
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supernovae type la. This is due to the significant non-Gaussianity of the likelihood func- 
tion with respect to We- In order to find the best-fit value of this parameter along with the 
best-fit values of the remaining cosmological parameters other datasets should be extensively 
analyzed. 

The goal of this paper is to derive constraints on the parameters of models with scalar 
fields using different current and near future data and to present observational constraints 
on cosmological models with classical and tachyonic scalar fields with barotropic equation of 
state as dark energy. 

The paper is organized as follows. In Section 2 we discuss the parametrization of 
barotropic EoS parameter by its current and early values, the evolution of the scale factor in 
single- and multicomponent models and the potentials of classical and tachyonic scalar fields 
with barotropic EoS. In Section 3 we present the observational constraints on the parameters 
defining the barotropic EoS obtained from the currently available data. In Section 4 we 
forecast the precision, with which the expected Planck data will be able to constrain the 
cosmological parameters of models with scalar field dark energy with barotropic EoS. The 
conclusion are given in Section 5. 



2 Scalar field models of dark energy with barotropic equation of state 

The background Universe is assumed to be spatially flat, homogeneous and isotropic with 
Friedmann- Robertson- Walker (FRW) metric of 4-space 

ds"^ = gijdx^dx^ = {ri){d'iif' — 6a/3dx°'dx^), 

where rj is the conformal time defined by dt = a{r])drj and a{rj) is the scale factor, normalized 
to 1 at the current epoch (here and below we put c = 1). The Latin indices i, j,... run from 
to 3 and the Greek ones are used for the spatial part of the metric: a, (3, .. = 1, 2, 3. 

We consider a multicomponent model of the Universe filled with non-relativistic particles 
(cold dark matter and baryons), relativistic particles (thermal electromagnetic radiation and 
massless neutrino) and minimally coupled dark energy. The dark energy is assumed to be a 
scalar field with either Klein-Gordon (classical, below: CSF) or Dirac-Born-Infeld (tachyonic, 
below: TSF) Lagrangian 

L,ias = X -U{<P), Ltach = -U{CWl-2X, (2.1) 

where U{(j)) and U (^) are the field potentials defining the model of the scalar field, X = 
(p-i(j)'^/2 and X = ^;j^'*/2 are kinetic terms. We assume the homogeneity of background scalar 
fields {(f){'x,i]) = 4>{'q), C(x, ?y) = S,{t])), so that their energy density and pressure depend only 
on time: 



Paas=X + U{(l)), Pdas = X-U{<P), (2.2) 
Ptach 



^J\-2X 

The EoS parameters tu^e = Pde/ Pde for these fields are 



^^^^ Ptach = -U{iWl-2X. (2.3) 



X — u 

Wclas = X + U ' ^^"""^ = 2X - I. (2.4) 
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Using the last relations the field variables and potentials can be presented in terms of densities 
and EoS parameters as: 



m-^o = ±j^ ^rwi^) ' (2-^) 

U{a) = P'^eia)[l-w{a)] ^^.6) 



for the classical Lagrangian and 



a 



da'y/l + w{a') 



Uia) = pdeia)y'-wia) (2.8) 

for the tachyonic case. 

The dynamics of expansion of the Universe is fully described by the Einstein equations 

R^, - Ur = SvtG (if) +T^;^+ T^'^) , (2.9) 



2^ -' \ V V 

where Rij is the Ricci tensor and t/'"^ Tj'^\ T^f'^'^ are the energy-momentum tensors of 

''J ''J ''J 

non-relativistic matter (m), relativistic matter (r), and dark energy (de) correspondingly. 
Assuming that the interaction between these components is only gravitational, each of them 
should satisfy the differential energy-momentum conservation law separately, which for a 
perfect fluid with density p„ and pressure pn related by the equation of state pn = WnPn 
yields: 

ap'^ = -?,Pn{l + Wn), (2.10) 

here and below a prime denotes the derivative with respect to the scale factor a. For the 
non-relativistic matter Wm = and pm = Pm^ , for the relativistic one Wr = 1/3 and 
Pr = Pr^^ CL~'^. Hereafter "0" denotes the current values. 

The EoS parameter w and the adiabatic sound speed = Pde/ Pde of dark energy are 
related by the differential equation: 

aw' = 'i{l + w){w - cl), (2.11) 

Note that in the case of studied types of dark energy the adiabatic sound speed is not the 
true velocity of sound propagation. The scalar fields have non-negligible entropy and thus 
non-adiabatic pressure perturbations. In the dark energy rest frame the total pressure per- 
turbation can be presented as 5pde = c^s^de-, where the effective sound speed = p^x / P-,x = 
L,x /{2XL,xx +L,x ) is different for the classical and tachyonic scalar fields (1 and —w cor- 
respondingly). Therefore, the adiabatic sound speed is a useful quantity having the form of 
corresponding thermodynamical quantity. Below we will see its real physical meaning. In 
general can be a function of time, but here we assume it to be constant: = const. In 
this case the time derivative of pde is proportional to the time derivative of pde or in integral 
form: 

Pde = clpde + C, (2.12) 
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where C is a constant. The above expression is the generahzed hnear barotropic equation of 
state. The solution of equation (2.11) for = const is 



(\ (l + cg)(l + ^o) ^ 

1 + Wq — [Wq — C^jfl'^*-^ 

where the integration constant wq is the current value of w. One can easily find that (2.13) 
gives (2.12) with C = p^^}{wo — c^), where p'^^ is the current density of dark energy. Substi- 
tuting (2.13) into (2.11) we see that for quintessence fields {wq > —1) the derivative of EoS 
parameter with respect to the scale factor is negative for > wq and positive for < wq. 

Therefore, we have two values wq and defining the EoS parameter w for any scale 
factor a. As it follows from (2.13), parameter corresponds to the EoS parameter at the 
beginning of expansion: We = w{0) = c^. Below we will call this parameter We instead of 
as such notation reflects better its physical meaning. The differential equation (2.10) with w 
in the form (2.13) has the analytic solution 

P^e = y YT^^ j ■ (2.14) 

Using the dependences of densities of each component on the scale factor the following equa- 
tions for background dynamics can be deduced from the Einstein equations (2.9): 



H = HQyJ^r/a^ + ^m/a? + ^d^f{a), (2.15) 
_ 1 2Q.r/a'^ + Vt^/a^ + (1 + ^w)Q.def{a) 



2 Qr/a^ + ^m/a^ + ^def{a) 



(2.16) 



where /(a) = [(\ + WQ)a~'^^^^^'''> +We — wq\/ {1 + We)- Here H = a/a? \s the Hubble parameter 
(expansion rate) and q = — (aa/a? — l) is the acceleration parameter (" ' "= d/drj). The 
equations (2.15)-(2.16) completely describe the dynamics of expansion of the homogeneous 
and isotropic Universe. 

In our previous paper [11] we have analyzed in detail three possible scenarios of the 
future evolution of the Universe. For We > wq the dark energy will tend to mimic a cosmo- 
logical constant in the future, such a Universe will expand forever as in de-Sitter inflation. 
For We = Wq (the simplest case) the future of the Universe is eternal power-law expansion. 
For We < Wq the dark energy turns away from its repulsive properties and in the future the 
expansion of the Universe will turn into contraction. 

For a realistic multicomponent model it is possible to find the time dependence of the 
scale factor from (2.15) only numerically. However let us consider the simple scalar field 
model of spatially-flat Universe fllled only with the scalar fleld with barotropic EoS. In this 
case the equation has the analytical solutions for the evolution of scale factor 

/ 1 + \ 3(l + »e) 2 / 3 , 

a{t) = sinh3(i+».) -^(l + 7x;g)(u;g-u'o)//oi 

\We-WoJ \2 

for We > Wq and 

1 

/ 1 + tfn \ 3{l+we) 2 / 3 , \ 

ait) = sin3(i+»e) -y/{l + We){wo-We)Hot 

\ Wo — We/ V 2 / 
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Figure 1. Evolution of the scale factor in cosmological models with scalar fields with barotropic 
EoS as dark energy. Black lines: We < wq, red: We > wq. Solid lines - numerical solutions for the 
multicomponent models, dotted - analytical solutions for the simple scalar field models. 

for We < wq. For We > wq this solution has been found previously in [14], the corresponding 
expression for We < wq can be easily obtained from it. 

In Fig. 1 the numerical solutions for the multicomponent model and the corresponding 
analytical ones for the simple scalar field model are presented for both cases: We > wq and 
We < Wq. We see that indeed there exist two possible scenarios for the future evolution 
of the Universe. In the case We > wq in far future the Universe will experience eternal 
asymptotically de-Sitter expansion, while in the case We < wq the cosmological expansion 
will slow down reaching the turnaround time after which the Universe will collapse. So, in 
the latter case the whole history of the Universe is limited in time. The difference between 
the corresponding curves in single-multicomponent models arises from the fact that in the 
multicomponent Universe at the early stages of expansion the relativistic and non-relativistic 
matter dominate. At scale factors corresponding to the dark energy domination in the 
multicomponent model the shapes of both curves become similar but are shifted in time. 

In the simple scalar field model from (2.5)-(2.8) it is easy to obtain the analytical 
expressions for the field potentials. For CSF they read: 



3Hq We - Wq 
SttG l + We 



+ 



SttG 1 + We 2 



l-WeWe-Wp 




in the case of We > wq and 
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Figure 2. The potentials of scalar fields with barotropic EoS. Blue lines: We < wq, black: We > wq. 
Solid lines - numerical solutions for multicomponent models, dotted - analytical solutions for simple 
scalar field models. Left: CSF, right: TSF. 



in the case We < wq. The potential of TSF with We > wq is: 



3ifo^ We 



Wo 



8ttG l + We 



sin ( -HQ^/we - Wo - ^o) + tan ^ 



sin^ ( -Ho^/we - Wo - Co) + tan ^ 



Wo 



l + wo 



We cos^ ( -HoVwe - Wo " Co) + tan ^ 



We 



Wo 



1 + Wo 



the corresponding potential of TSF with We < wq reads: 



t^(e-Co) 



sinh 



-We COSh^ 



SH^ wo 



We 



inG 1 + We 



sinh ( -HoVwq - We - Co) 



tanh 



-Ho^/w^^^^e - Co) + tanh ^ 



-Ho\fwo - We (C - Co) + tanh ^ 



Wo 



We 



l + wo 

Wq - We 
l + Wo 



We 



Wo 



l + Wo 



Wo 



We 



l + Wo 



For We > Wo these potentials have been found in [14], for We < wq the corresponding poten- 
tials can be derived from them. In Fig. 2 the potentials of CSF and TSF are presented for the 
realistic multicomponent model as well as for the simple scalar field one. The potentials in 
both models are very different at early epoch but become similar in the epoch corresponding 
to the dark energy domination in the multicomponent model. In both models the potentials 
of CSF with We < Wo become negative at some time in future while the potentials of TSF 
become imaginary. The potentials of fields with We > wq have no such peculiarities. It should 
be noted that for scalar fields with We > wo the field variables tend to a finite value at infinite 
time. 



3 Observational constraints from current datasets 

In the previous section it has been shown that the dynamics of the expansion of the Universe 
and its future are determined by the relation between the parameters wq and We- To find 
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Figure 3. From top to bottom: the relative differences of CMB temperature fluctuations power 
spectra, matter density power spectra and distance moduli in cosmological models with CSF with 
barotropic EoS as dark energy and in the fiducial ACDM model. Black lines correspond to wq = —0.99, 
blue to Wo = —0.77, green to wq = —0.55 and red to wq = —0.33. Solid lines correspond to We = —0.99, 
dashed ones to We = —0.66, dash-dotted ones to We = —0.33 and dash-3-dotted ones to We = 0. 



out which scenario is valid these parameters should be determined from the observational 
data. As the values of other cosmological parameters are unknown, the determination has to 
be performed for the full set of cosmological parameters which involves also the dark energy 
density parameter i^dej the physical density parameters of baryons Oft/i^ and cold dark matter 
^cdmh'^, the Hubble constant Hq, the spectral index of initial matter density power spectrum 
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Figure 4. One-dimensional marginalized posteriors (solid lines) and mean likelihoods (dotted 
lines) for the main cosmological parameters derived from the datasets SN Union2, SN SDSS 
SALT2, SN SDSS MLCS2k2, SN Union2+HST+BBN, SN SDSS SALT2+HST+BBN and SN SDSS 
MLCS2k2+HST+BBN (from top to bottom). Black lines - ACDM model, red - model with CSF. 



Ug, the amplitude of initial matter density power spectrum Ag and the reionization optical 
depth r (here and below h = Hq/WO km/[s ■ Mpc]). Instead of Hq we vary the commonly 
used parameter 6 defining the ratio of the sound horizon to the angular diameter distance 
multiplied by 100. These are nine unknown parameters, but the number of independent ones 
is 8, since we assume spatial flatness, hence the dark energy density parameter is obtained 
from the zero curvature condition: fi^e = ^ — — ^cdm- We neglect the contribution from 
the tensor mode of perturbations and the masses of active neutrinos. 

3.1 Method and data 

In our previous paper [11] we have found that all cosmological parameters are determined well 
with the exception of We, which remains essentially unconstrained by the combined datasets 
used there. In order to find the best-fit value of this parameter and its confidence limits 
for scalar fields with both types of Lagrangians in this paper we perform the Markov Chain 
Monte-Carlo (MCMC) analysis for different combined datasets. 
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For our analysis we use the publicly available package CosmoMC [15, 16] including the 
publicly available code CAMB [17, 18] for the calculation of the model predictions. This 
code has been modified to include the dark energy models proposed here as described in our 
previous paper [11]. 

Each of the performed MCMC runs has 8 chains converged to i? — 1 < 0.01. 

We use the following datasets: 

• CMB temperature fluctuations and polarization angular power spectra from the 7-year 
WMAP observations (hereafter WMAP7) [19, 20]; 

• Baryon acoustic oscillations in the space distribution of galaxies from SDSS DR7 (here- 
after BAO) [21]; 

• Power spectrum of luminous red galaxies from SDSS DR7 (hereafter SDSS LRG7) [22] - 
in this case we obtain the nonlinear correction of the small-scale matter power spectrum 
using the version of halof it modified to include the background dynamics in models 
with scalar fields with barotropic EoS neglecting the dark energy perturbations which 
have been found to be significantly smaller than the dark matter ones at these scales 
at the current epoch [11, 23]; 

• Hubble constant measurements from HST (hereafter HST) [24]; 

• Big Bang Nucleosynthesis prior on baryon abundance (hereafter BBN) [25, 26]; 

• Union2 supernovae la compilation including 557 SN with SALT2 method of light curve 
fitting (hereafter SN Union2) [27]; 

• SDSS supernovae la compilation (hereafter SN SDSS) [28] - full sample includes 288 
SN; both SALT2 [29] and MLCS2k2 [30] (modified by [28]) methods of light curve 
fitting are used. 

3.2 The effect of scalar fields vi^ith barotropic EoS on characteristics of expansion 
dynamics and large scale structure of the Universe 

Let us briefly discuss the effect of the dark energy models studied here on measurable 
characteristics of the dynamics of expansion and large scale structure of the Universe. In 
Fig. 3 we present the relative differences of the CMB temperature fluctuations power spectra 
(C7^-Cf ^^*^V<^f ^^^'^^ matter density power spectra {P{k)- P^^^'^^^\k))/p(-^^^^^\k) 
and distance moduli {fi — ^('^'-^^^^))//i(^'^^*'^) (the distance modulus is given by the lumi- 
nosity distance as /x = 51og]^o + 25) for models with CSF with different values of wq and 
We- For comparison we use the fiducial flat ACDM model with parameters Ofo/i^ = 0.0223, 
^cdmh"^ = 0.105, h = 0.7, r = 0.09, = 0.95 and As = 2.0 • 10"^. The main cosmological 
parameters in the models with a scalar field are taken to be equal to the fiducial ones. The 
symbols in the figure show the relative errors of the used observational data. 

In the top panel we see that the dark energy with wq = —0.99 leads to differences in 
CMB temperature fluctuations power spectrum which are at the percent level. On the other 
hand the dark energy with We = causes at low multipoles relative differences of factors 
~ 1.5 — 4 for all values of wq which are not very close to —1. This suggests that the values 
of We close to should be ruled out by the CMB data with sufhciently high precision. The 
relative errors of observational WMAP7 data points are smaller than most of the relative 
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differences between model spectra up to £ ~ 1000, however tlie cosmological parameters, 
wliich have been fixed, in Fig. 3 should in practice be determined from the data jointly with 
the dark energy parameters wq and We ■ This reduces the possibility of distinguishing between 
the models with different values of wq and We ■ 

The CMB allows a reliable determination of the main cosmological parameters: we have 
checked that the one-dimensional marginalized posteriors and mean likelihoods derived from 
WMAP7 only for the ACDM model and the model with CSF with barotropic EoS are close, 
which means that the proposed dark energy model does not cause the large deviations of the 
best-fit values of main cosmological parameters from their values in models with cosmological 
constant, however for We the shapes of the marginalized posterior and mean likelihood are 
different, which reflects the fact that this parameter cannot be constrained using CMB data 
only due to the small dark energy density (comparing to the critical one) at recombination. 

From the middle panel of Fig. 3 it follows that the dark energy with barotropic EoS also 
causes differences in ^(A;) at level of percents and tens of percents for all values of wq and 
We not very close to —1. This suggests that the data on the matter power spectrum could be 
useful for constraining the parameters of the studied models. Note that in the case of P{k) 
the maximal relative differences are smaller than for the CMB and their absolute value never 
exceeded 1. Similarly to the CMB spectrum from WMAP7, for values of wq not very close 
to —1 the relative errors of the SDSS LRG7 data are smaller than the relative differences of 
power spectra between models with CSF with barotropic EoS and ACDM. 

The data on the power spectrum of luminous red galaxies should be used combined with 
CMB datasets or priors, in such combinations they serve as an additional constraints on the 
parameters of the scalar fields with barotropic EoS and other cosmological parameters. 

In the bottom panel the relative differences of the distance moduli are shown. They 
decrease with redshift while the distance modulus and luminosity distance grow. These 
differences are the smallest among all studied, they do not exceed 0.03, but up to z ~ 0.5 
they are larger than most of the relative errors of supernovae distance moduli from used 
compilations and up to z ~ 1 they remain comparable to these relative errors of data for all 
values of the current scalar field EoS parameter that are not very close to —1. The supernovae 
data are known to be the key data for determination of the best-fit values of cosmological 
parameters and their confidence range in models with dark energy with variable EoS. At 
the beginning of the next subsection we shall check their importance for constraining our 
particular model. 

Other data providing distance information are also BAO and HST. BAO should be 
used in combination with WMAP7 or with additional priors on parameters derived from the 
CMB datasets, in these combinations they can provide additional constraints on values of 
all cosmological parameters. HST data are useful for the determination of Hq as well as to 
provide an additional constraint on deviations of cosmological models from the fiat ACDM 
model with Q\ = 0.7 and Hq = 74.2 km/[s ■ Mpc] at z = 0.04 when combined with other 
datasets. We also use the BBN prior on baryon abundance that puts an additional constraint 
on Qiyh? . 

In this paper we use fiat priors for all parameters, for wq and We the allowed ranges are set 
to —1 <wq < —0.33 and —l<'We<0. The reason for the lower bound on both parameters is 
that we study the quintessential scalar fields, for which the field variables become imaginary 
if the EoS parameter is smaller than —1. The upper bound on wq corresponds to the value 
at which the accelerated expansion of the Universe at current epoch is no longer possible. 
The upper bound on the early EoS parameter can be obtained using e.g. Fig. 3, in the top 
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Figure 5. Onc-dimcnsional marginalized posteriors (solid lines) and mean likelihoods (dotted lines) 
for the datasets WMAP7+HST+BBN, WMAP7+HST+BBN+BA0, WMAP7+HST+BBN+SDSS 
LRG7, WMAP7+HST+BBN+SN Union2 and WMAP7+HST+BBN+SN SDSS SALT2 (from top to 
bottom). From left to right: the functions for the dark energy density ^de-, the current value of dark 
energy EoS parameter wa for both CSF (black lines) and TSF (red), the early value of dark energy 
EoS parameter We for CSF and the early value of dark energy EoS parameter We for CSF. 

panel of which we see that We = can cause so large changes in the CMB power spectrum 
that larger values of We should be avoided. 

3.3 Results 

At the beginning we analyse the possibility of constraining the cosmological parameters of 
models with scalar fields with barotropic EoS as dark energy using the supernovae data. 

In 3 upper rows of Fig. 4 we present the results of determination of the main cosmo- 
logical parameters along with wq and We from the supernovae data only. We see that for 
the parameters defining the barotropic EoS, wq and We, the one-dimensional marginalized 
posteriors and mean likelihoods are relatively close. This means that the supernovae data 
indeed are important for constraining the studied dark energy models. However, the param- 
eters VLiyh?, Qcdmh^ Si,nd 9 remain essentially unconstrained for both ACDM model and model 
with barotropic EoS. In 3 lower rows of Fig. 4 we present the one-dimensional marginalized 
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Figure 6. Two-dimensional mean likelihood distributions in the plane We — wo for the same datasets 
and models as in Fig. 5. Solid lines show the la and 2a confidence contom'S. 

posteriors and mean likelihoods for the same models and supernovae datasets combined with 
BBN and HST. We see that now the parameters of ACDM model are constrained well, while 
for the model with scalar field ^l^h'^ is constrained as well as in the ACDM case but Qcdm.h'^ 
and 6 are constrained better than by supernovae data only though not perfectly. All 3 su- 
pernovae datasets provide the constraints with comparable level of accuracy, however it is 
worth to note that SN SDSS with MLCS2k2 fitting puts the clearly weakest constraint on 

Wo- 

Now we analyse the constraints on studied dark energy models from the combined 
datasets including WMAP7 data. In the first step we analyze 6 combined datasets: 



- 12 - 



• WMAP7+HST+BBN, 



• WMAP7+HST+BBN+BA0, 

• WMAP7+HST+BBN+SDSS LRG7, 

• WMAP7+HST+BBN+SN Union2, 

• WMAP7+HST+BBN+SN SDSS SALT2 and 

• WMAP7+HST+BBN+SN SDSS MLCS2k2. 

In Fig. 5 we see that for the datasets WMAP7+HST+BBN, WMAP7+HST+BBN 
+BAO, WMAP7+HST+BBN+SDSS LRG7, WMAP7+HST+BBN+SN Union2 and 
WMAP7+HST+BBN+SN SDSS SALT2 the shapes of the marginalized posterior distri- 
butions and the mean likeHhoods for the early value of EoS parameter are different, this 
indicates a significant non-Gaussianity of the likelihood function for Wg. The difference be- 
tween two-dimensional marginalized posteriors and mean likelihoods in the plane We — wq, 
which is shown in Fig. 6, confirms this conclusion. As it has been shown in our paper [13] this 
is also the case for the combination WMAP7+HST+BBN+SN Union with the light curves 
of supernovae from the Union compilation [31] fitted using the SALT method [32]. 

The best-fit values of the cosmological parameters (obtained from the best-fit sample) 
and their la limits from the extremal values of the N-dimensional distribution are pre- 
sented in Tables 1-2 for the combined datasets WMAP7+HST+BBN, WMAP7+HST+BBN 
+BAO, WMAP7+HST+BBN+SDSS LRG7, WMAP7+HST+BBN+SN Union2 and 
WMAP7+HST+BBN+SN SDSS SALT2. 

The situation is different for the dataset WMAP7+HST+BBN+SN SDSS MLCS2k2. 
In this it can be seen in the top panels of Fig. 7, the shapes of mean likelihoods 

and the one-dimensional marginalized posteriors for We are similar to the half-Gaussian with 
center at the boundary of the allowed range of values, We = —1. This dataset can be used 
for a reliable estimation of the early EoS parameter. The difference between both curves is 
a signal of non-Gaussianity, which is however substantially reduced compared to all other 
datasets. The two-dimensional mean likelihoods and marginalized posteriors in the plane 
We — Wq presented in the top panels of Fig. 8 support this conclusion. The shapes of the 
high-likelihood regions are similar to the shapes of la and 2a confidence contours. 

Hence, we have found that the data on SN la from the SDSS compilation with modified 
MLCS2k2 fitting of light curves allow to constrain We while the same data with SALT2 fitting 
do not. This is a demonstration of the well-known discrepancy between SALT2 and MLCS2k2 
which is due mainly to the different rest-frame U-band models and the assumptions about 
the color variations in both fitting methods [28]. 

In paper [13] we have performed similar MCMC runs for the combined datasets in- 
cluding the SN subset NEARBY+SDSS (136 SN) from the SDSS compilation, for which this 
discrepancy is smallest [28] . In this case the parameter We remains unconstrained for SN data 
with both light curve fitting methods. Therefore, the non-Gaussianity of the likelihood func- 
tion with respect to We is reduced by inclusion of the higher-redshift SN samples, for which 
the treatment of the MLGS2k2 method differs significantly from the corresponding treatment 
of the SALT2 method. In Fig. 9 we present the one- and two-dimensional marginalized pos- 
teriors and mean likelihoods derived from the MCMG runs performed for the case of GSF 
for combined datasets including WMAP7, HST, BBN, NEARBY+SDSS SN data and in 
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Table 1. The best-fit values for cosmological parameters and the la hmits from the extremal values of 
the N-dimensional distribution determined for the case of CSF by the MCMC technique from the com- 
bined datasets WMAP7+HST-hBBN, WMAP7+HST-FBBN+BA0, WMAP7+HST+BBN+SDSS 
LRG7, WMAP7-hHST+BBN-hSN Union2 and WMAP7+HST-hBBN+SN SDSS SALT2. 

addition different higher-redshift subsamples of the SDSS compilation with MLCS2k2 light 
curve fitting. We see that the reduction of the non-Gaussianity is due mainly to the inclusion 
of the SNLS subsample, for which the difference between SALT2 and MLCS2k2 is found to 
increase with redshift [28]. 

Let us fix the parameters Q^h'^, Qcdm^'^ ^ to their fiducial values and determine the 
constraints on wo and We from the supernovae data only. In the 2 upper rows of Fig. 10 the 
one-dimensional marginalized posteriors and mean likelihoods are shown for We and wo- We 
see that in this case none of the supernovae data allow a reliable determination of We- The 
constraints on wo obtained from SN SDSS MLCS2k2 are also weak. Now we combine the 
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Table 2. The best-fit values for cosmological parameters and the la hmits from the extremal values of 
the N-dimensional distribution determined for the case of TSF by the MCMC technique from the com- 
bined datasets WMAP7+HST-hBBN, WMAP7+HST+BBN+BA0, WMAP7+HST+BBN+SDSS 
LRG7, WMAP7-hHST+BBN-hSN Union2 and WMAP7+HST-hBBN+SN SDSS SALT2. 

supernovae data with the HST ones. In the 2 bottom rows of Fig. 10 the corresponding the 
one-dimensional marginalized posteriors and mean likelihoods for We and wo are presented. 
We find that the SN SDSS SALT2 and SN Union2 data still do not allow to put constraints 
on We while in the case of SN SDSS data the marginalized posterior and mean likelihood are 
close. Thus the SN SDSS MLCS2k2 data constrain the values of We when combined with 
the additional constraint on the deviations of model from ACDM coming from HST. The 
transition of positions of the peaks of posterior and likelihood from values close to in Fig. 4 
to —1 is due to the fact that the parameters Qcdmh'^ and 6 which are highly correlated with 
We and wq in the case presented in Fig. 4 are now fixed to their fiducial values. The current 
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Figure 7. Onc-dimensional marginalized posteriors (solid lines) and mean likelihoods (dotted ones) 
for the combined datasets WMAP7+HST+BBN+SN SDSS MLCS2k2, WMAP7+HST+BBN+SN 
SDSS MLCS2k2+BAO and WMAP7+HST+BBN+SN SDSS MLCS2k2+SDSS LRG7 (from top to 
bottom). From left to right: the functions for the dark energy density fide, the current value of dark 
energy EoS parameter wq for both CSF (black lines) and TSF (red), the early value of dark energy 
EoS parameter We for CSF and the early value of dark energy EoS parameter We for CSF. 

value of the dark energy EoS parameter wq is constrained well by all 3 supernovae datasets 
combined with HST. While SN SDSS SALT2 and SN Union2 prefer values close to —1, the 
best-fit values of wq obtained from SN SDSS MLCS2k2 are closer to -0.8. 

A general discussion of the differences, benefits and limitations of the SALT2 and 
MLCS2k2 light curve fitting methods is beyond the scope of this paper. Here we should 
mention only that the MLCS2k2 is believed to be more model-independent than SALT2 
which fits the global parameters M, a and /3 determining the distance moduli together with 
the cosmological ones. Let us compare the root mean squares of errors of different super- 
novae samples. For SN Union2 it is 0.269, for SN SDSS SALT2 it equals 0.260 while for 
SN SDSS MLCS2k2 0.233. Thus the supernovae from the SDSS compilation with MLCS2k2 
fitting have the smallest errors among all used compilations. These small difference results in 
the possibility of reliable determination of the best-fit value of We and its confidence ranges. 
Considering the subsamples of the SN SDSS compilation, the datasets NEARBY+SDSS SN 
have the root mean square of errors 0.198 for MLCS2k2 and 0.191 for SALT2 - these values 
are almost the same, so there is no difference in possibility of constraining Wg from these data. 
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Figure 8. Two-dimensional mean likeliliood distributions in the plane We — wq for corresponding 
datasets and models from Fig. 7. Solid lines show the la and 2a confidence contours. 

Among the higher-redshift subsamples of SN SDSS compilation the samples ESSENCE, HST 
and SNLS have the root mean squares of errors 0.261, 0.283 and 0.246 correspondingly. Thus, 
the inclusion of SNLS subsample is crucial for constraining of the early value of dark energy 
EoS parameter since it has the smallest errors among all higher-redshift subsamples of SN 
SDSS compilation with MLCS2k2 light curve fitting method. 

As the next step we perform the similar MCMC runs for two combined datasets including 
WMAP7+HST+BBN+SN SDSS MLCS2k2 and the data on large-scale structure: 

• WMAP7+HST-FBBN+SN SDSS MLCS2k2-hBAO and 

• WMAP7+HST+BBN+SN SDSS MLCS2k2+SDSS LRG7. 

The one- and two-dimensional marginalized posteriors and mean likelihoods for these datasets 
are shown in middle and bottom panels of Fig. 7-8. These datasets allow the constraints on 
the value of We at comparable level of accuracy as the previous dataset WMAP7-|-HST-|-BBN 
+SN SDSS MLCS2k2. For the dataset WMAP7+HST+BBN+SN SDSS MLCS2k2+BAO 
the one-dimensional marginalized posterior and mean likelihood have the shape of the half- 
Gaussian with center at the boundary of the allowed range of values {we = — 1). The small 
difference between both curves indicates the slight non-Gaussianity of the likelihood function 
with respect to We, which however does not reduce the possibility of a reliable estimation of 
We from these data. For the set WMAP7-hHST-hBBN-FSN SDSS MLCS2k2-hSDSS LRG7 
they have the shape of the Gaussians with centers at slightly different and larger values of 
We- This difference signals that some non-Gaussianity of the likelihood function with respect 
to the early EoS parameter exists, however, as it is relatively small, we conclude that the 
last dataset can also be used for the reliable estimation of We- 

The two-dimensional We — wq marginalized posteriors and mean likelihoods presented 
in the middle and bottom panels of Fig. 8 support these conclusions. For the dataset 
WMAP7-FHST-FBBN+SN SDSS MLCS2k2-FBAO the shapes of the high-likelihood regions 



and of the la and 2a confidence contours are similar. For the combination WMAP7+HST 
+BBN+SN SDSS MLCS2k2+SDSS LRG7 the high-likelihood region lays partially outside 
the 2a confidence contour, this is a signal of the above mentioned non-Gaussianity, which 
can be however neglected since it is relatively small. 

From Fig. 7 we see that the fields mimicking a cosmological constant are excluded 
at the la confidence level for all 3 datasets for both Lagrangians. Moreover, for the set 
WMAP7+HST+BBN+SN SDSS MLCS2k2+SDSS LRG7 both such fields lay even shghtly 
outside the 2a confidence contour. This is due to the inclusion of SN SDSS data with light 
curves fitted using the modified MLCS2k2 method. The values of We close to are excluded 
nearly at the 2a confidence level for both fields and all 3 datasets. 

The best-fit values of the cosmological parameters and their la limits from the extremal 
values of the N-dimensional distribution are presented in Table 3 for CSF and TSF models 
with barotropic EoS for the combined datasets WMAP7+HST+BBN+SN SDSS MLCS2k2, 
WMAP7+HST+BBN+SN SDSS MLCS2k2+BAO and WMAP7+HST+BBN 
+SN SDSS MLCS2k2+SDSS LRG7. Note that these limits are significantly wider than 
the corresponding limits obtained from the one- and two-dimensional marginalized distribu- 
tions. We see that for all cases including SN SDSS MLCS2k2 the best-fit model has Wg < wq, 
thus the repulsive character of the scalar fields will stop and the expansion of such Universe 
will turn into collapse. 

It should be noted that the fields with classical and tachyonic Lagrangians cannot be 
distinguished by the currently available data: the differences between the best-fit parameters 
are within the corresponding la confidence limits (see also [12]). Moreover, as we can see 
in the bottom rows of Tables 1-3, the values of — log L differ by less than 1 for the best-fit 
parameters of models with CSF and TSF obtained from the same datasets. Note that the 
values of -logL for the WMAP7+HST+BBN+SN SDSS SALT2 dataset are higher than 
those for the same dataset including SN SDSS with MLCS2k2 light curve fitting method. 

It should be noted that, as it can be seen in Fig. 7 and 8, the dark energy density Qde 
and its current value of EoS parameter wq are well constrained by all the used combined 
datasets. The curves for CSF and TSF overlap. 

Finally, let us discuss the best-fit values of Hubble constant obtained from different 
datasets. As it can be seen in Tables 1-3, for the combined datasets, which do not include 
SN SDSS MLCS2k2 data, the best-fit values of Hq are in the range 68.2-72.3 km/{s ■ Mpc) 
which is closer to Hq = 74.2 km/{s ■ Mpc) from [24] than the range 65.9-67.1 km/{s ■ 
Mpc) obtained when we include these data. Performing the MCMC runs for the model 
with CSF and the combined datasets including WMAP7, HST, BBN and different subsam- 
ples of SN SDSS compilation with MLCS2k2 light curve fitting we have found the follow- 
ing fit values of Hubble constant (in km/{s ■ Mpc)): 71.8 for NEARBY+SDSS SN (and 
71.5 for the SALT2 fitting method), 70.1 for NEARBY+SDSS+ESSENCE SN, 68.5 for 
NEARBY+SDSS+HST SN, 67.7 for NEARBY+SDSS+SNLS SN and 67.5 for NEARBY 
+SDSS+ESSENCE+HST SN. Thus, in the case of MLCS2k2 fitting of light curves the best- 
fit value of Hq is lowered mainly by the subsamples either ESSENCE+HST or SNLS from 
the SN SDSS compilation. 

4 Error forecasts for the Planck experiment 

In the previous section we have determined the observational constraints on cosmological 
parameters in models with scalar fields with barotropic EoS. Now we are going to discuss 
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Figure 9. Top: one-dimensional marginalized posteriors (solid lines) and mean likelihoods (dotted 
ones) for CSF and the combined datasets WMAP7+HST+BBN+SN SDSS MLCS2k2 for the sub- 
samples NEARBY+SDSS+ESSENCE, NEARBY+SDSS+HST, NEARBY+SDSS+ESSENCE-I-HST 
and NEARBY-hSDSS+SNLS of SN SDSS compilation (from left to right). Bottom: corresponding 
two-dimensional mean likelihood distributions in the plane We — wq. Solid lines show the la and 2a 
confidence contours. 

the precision, with which the expected Planck data on CMB anisotropics will allow us to 
estimate the parameters determining the barotropic EoS. 

If the shape of the likelihood function cannot not be safely assumed to be Gaussian, the 
most reliable forecasting technique is the full MCMC analysis of mock data. To generate a 
Planck mock dataset we have used the publicly available code FuturCMB [33] . The method of 
mock data generation and all necessary modifications of CosmoMC are thoroughly described 
in [34]. The fiducial Cgs have been computed by CAMB for the set of best-fit parameters 
obtained for CSF from the dataset WMAP7-hHST-hBBN-FSN SDSS MLCS2k2-FSDSS LRG7 
(Table 3, column 3). The seed has been chosen as 150. The generated mock dataset involves 
the CMB temperature fluctuations, polarization and the weak lensing deflection angle power 
spectra. 

We assume that the Planck experiment has 3 channels, for which we choose in turn 
Ofwhm, (Tt and aE to be 9.5 arcmin, 6.8 per pixel and 10.9 jiK per pixel; 7.1 arcmin, 6.0 
fiK per pixel and 11.4 nK per pixel; 5.0 arcmin, 13.1 per pixel and 26.7 per pixel 
correspondingly. The observed sky fraction is taken to be fsky = 0.65. 

We have performed MCMC runs similar to those in the previous section and found that 
the Planck data alone as well as the combinations Planck -|-HST-|-BBN, Planck -|-HST-|-BBN 
+BAO, Planck-hHST-hBBN-hSDSS LRG7, Planck-hHST-hBBN-FSN Union2 and Planck 
-|-HST-|-BBN-|-SN SDSS SALT2 do not allow a reliable determination of the early value 
of EoS parameter. The likelihood for We is significantly non-Gaussian for any precision of 
the CMB data when the SN SDSS data with MLCS2k2 light curve fitting are not included. 

In Fig. 11-12 the one- and two-dimensional marginalized posteriors and mean likelihoods 
are presented for the datasets Planck-^BBN-^HST-FSN SDSS MLCS2k2, Planck+BBN-hHST 
-hSN SDSS MLCS2k2-FBAO and Planck-^BBN-^HST-FSN SDSS MLCS2k2-hSDSS LRG7. 
Also in this case, the higher-redshift SN from the full SDSS compilation with the light curves 
fitted by MLCS2k2 reduce the non-Gaussianity of the likelihood for We- 
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Figure 10. Onc-dimcnsional marginalized posteriors (solid lines) and mean likelihoods (dotted lines) 
for the datasets SN SDSS MLCS2k2, SN SDSS SALT2, SN Union2 (2 upper rows, from left to right) 
and SN SDSS MLCS2k2+HST+BBN, SN SDSS SALT2+HST+BBN, SN Union2+HST+BBN (2 
lower rows, from left to right) for the models where only wq and We are allowed to vary. 

The substitution of the 7-year WMAP data on CMB anisotropy by the Planck mock 
dataset reduces significantly the errors, so that the values of We close to appear to be 
far beyond the 2a marginalized confidence contours. It is worth noting that from the two- 
dimensional marginalized distributions it follows that the A-term dark energy can be excluded 
at the 2a confidence level for both Lagrangians and all datasets. 

The best-fit values of the cosmological parameters for models with both fields and the 
corresponding la limits from the extremal values of the N-dimensional distribution are pre- 
sented in Table 4. The best-fit values are close to the values obtained from the corresponding 
datasets including WMAP7 data within la limits, as expected. All best-fit models have 
We < Wq as the fiducial model. Note that for the parameters Q^h^, Qcdm^'^^ n-s, log (10^*^^^), 
Zrei and to the presented in Table 4 la uncertainties are few times smaller than the corre- 
sponding uncertainties presented in Table 3. The uncertainties of determination of We are 
significantly reduced by inclusion of the Planck mock data. The corresponding uncertainties 
of determination of the other dark energy parameters, il.de and wq, as well as of the Hubble 
constant Hq are also smaller than presented in Table 3 ones. 

Finally we want to check the reliability of the forecast. For this purpose we generate 4 
additional independent Planck mock datasets with different seeds: 50, 100, 200 and 250. We 
have performed MCMC runs for these additional datasets combined with SDSS LRG7, SN 
SDSS MLCS2k2, HST and BBN. In Fig. 13 the best-fit values and la fimits from the extremal 
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Table 3. The best- fit values for cosmological parameters and the \a hmits from the extremal values 
of the N-dimensional distribution determined by the MCMC technique from the combined datasets 
including SN SDSS data with light curve fitting MLCS2k2 as well as HST and BBN. The current 
Hubble parameter Hq is in units km/{s ■ Mpc), the age of the Universe to is given in Giga years. 



values of the N-dimensional distribution are shown for the parameters wq and We- We see 
that the best-fit values obtained from all datasets are within the la confidence limits and the 
limits are generally consistent with each other. It can be stated with high confidence that 
the values We > —0.1 should be excluded by the combined datasets including forthcoming 
Planck data. This is consistent with our conclusion that the models with values of Wg close 
to could possibly be distinguishable from the corresponding models with We close to —1 by 
the Planck data [12, 13]. Note that from Fig. 11-12 and Table 4 it can be deduced that CSF 
and TSF cannot be distinguished by CMB data from the next generation experiments since 
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Figure 11. One-dimensional marginalized posteriors (solid lines) and mean likelihoods (dotted lines) 
for the combined datasets Planck+HST+BBN+SN SDSS MLCS2k2, Planck+HST+BBN+SN SDSS 
MLCS2k2+BAO and Planck+HST+BBN+SN SDSS MLCS2k2+SDSS LRG7 (from top to bottom). 
From left to right: the functions for the dark energy density fide, the current value of dark energy 
EoS parameter wq for both CSF (black lines) and TSF (red), the early value of dark energy EoS 
parameter We for CSF and the early value of dark energy EoS parameter We for CSF. 

the differences between the models with both fields are smaller than the la confidence limits. 
In the bottom row of Tables 4 we see that the values of — log L differ for less than 1 for the 
best-fit parameters of models with CSF and TSF obtained from the same datasets. As it 
can be seen in Fig. 11, the dark energy density fide and its current value of EoS parameter 
Wo are constrained well by all used combined datasets including Planck mock data and the 
curves for CSF and TSF overlap. 

5 Conclusion 

We have constrained the parameters of cosmological models with classical and tachyonic 
scalar fields with barotropic equation of state as dark energy using combined datasets in- 
cluding the CMB power spectra from WMAP7, the Hubble constant measurements, the 
Big Bang nucleosynthesis prior, the baryon acoustic oscillations in the space distribution of 
galaxies from SDSS DR7, the power spectrum of luminous red galaxies from SDSS DR7 and 
the light curves of SN la from 2 different compilations: Union2 (SALT2 light curve fitting) 
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Figure 12. Two-dimensional mean likeliliood distributions in the plane We — wq for corresponding 
datasets and models from Fig. 11. Solid lines show the la and 2a confidence contours. 
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Figure 13. The best-fit values of We and wq along with la limits from the extremal values of the 
N-dimensional distribution for 5 independent Planck mock datasets (with seeds 50, 100, 150, 200 and 
250 from left to right in each panel). 



and SDSS (SALT2 and MLCS2k2 light curve fittings). We have found that the parameter 
corresponding to the value of w at early times, is essentially unconstrained by most of the 
currently available data due to the significant non-Gaussianity of the likelihood function for 
We- To determine the best-fit value and the la confidence ranges of We the combined datasets 
including SN data from the full SDSS compilation with MLCS2k2 fitting of light curves have 
to be used, since only these SN data reduce the non-Gaussianity sufficiently. In these cases 
the best-fit scalar fields have increasing EoS parameters, their repulsion properties recede 
and the Universe turns into contraction. 

We have also forecasted the uncertainties of the estimation of cosmological parameters 
of the studied models from the combined datasets including the data from the Planck ex- 
periment. We were especially interested in the precision, with which the Planck data will 
constrain the early EoS parameter value. We have found that the non-Gaussianity of the 
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Table 4. The best-fit values of cosmological parameters and Icr limits from the extremal values 
of the N-dimcnsional distribution determined by the MCMC technique from the combined datasets 
including SN SDSS data with light curve fitting MLCS2k2 and Planck mock data instead of WMAP7. 
All datasets include also HST and BBN. 



likelihood function with respect to w^ is not reduced by the expected Planck data alone. 
For the combined datasets including Planck mock data and SN data from the full SDSS 
compilation with MLCS2k2 light curve fitting method it is concluded that the models with 
We > —0.1 should be excluded at the 2a confidence level. 
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